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Abstract 

The computational Projective Dynamics method is used to analyze simulations of 
magnetization reversal in nanoscale magnetic pillars. It is shown that this method 
can be used to determine the magnetizations corresponding to the metastable min- 
imum and saddle point in the free energy, and the free-energy barrier associated 
with those points. For the nanopillars studied here, entropy is found to provide 
a significant contribution to the free-energy barrier which determines the reversal 
time scale. 
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Understanding magnetization reversal is important for technological applica- 
tions where the magnetic orientation of nanoscale regions must be quickly 
assigned. For information storage, two configurations of the magnetization are 
used to encode the state of a bit of data. During the assignment process, strong 
fields are applied to create a free-energy minimum for only one of these config- 
urations, and the bit almost certainly assumes the equilibrium configuration. 
For weaker fields, the two configurations correspond to local free-energy min- 
ima separated by a free-energy maximum. The minimum corresponding to a 
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magnetization parallel to the applied field is truly stable, while the antipar- 
allel minimum is higher in energy and therefore metastable. In this situation 
the configuration depends not only on the relative weights of the minima, 
but for device and human time scales, it also depends on the history of the 
configuration. At some point along the most probable path between the min- 
ima there must occur a free-energy maximum in the form of a saddle point. 
The free-energy difference between the maximum and either minimum, often 
called a free-energy barrier, determines the time scale for transitions of the 
configuration from that minimum to the other. 

To a first approximation these barriers and the curvature of the free energy 
near the extrema control the nonequilibrium dynamics of a system. Thus lo- 
cating free-energy extrema is essential for developing a detailed understanding 
of the dynamics of metastable states, and that understanding is important for 
such things as maximizing data integrity and enabling the technology of hybrid 
recording [1] , which uses lower-than-coercive applied fields to assign magnetic 
orientations in high-coercivity magnetic materials. While the free-energy min- 
ima can be easily located, the saddle point has proven much harder to measure 
[2]. Here we present results using the Projective Dynamics method [3,4,5] to 
probe the magnetization reversal of high- aspect-ratio nanoscale model mag- 
nets. 

It has already been shown that the Projective Dynamics method can be used 
to locate the saddle point [4,6,7]. The method involves projecting the original 
description of the dynamics in terms of a large number of variables onto a 
stochastic description in terms of one variable. For example, the dynamics of 
the thousands of individual spins in a nanoscale pillar can be projected onto 
the stochastic dynamics of the total magnetization along the long axis of the 
nanomagnet, M z . The transition rates between values of M z are measured 
by the probabilities P gr0 w, which correspond to an increase in the volume 
of stable magnetization, and P s hrmk, which correspond to a decrease in the 
volume. For values of M z with -P gr ow>-Pshrink, on average the volume of the 
stable magnetization grows. This corresponds to a negative local slope for 
the free energy. Likewise, -P gr0 w<-Pshrink corresponds to a positive local slope. 
Values of M z for which P gr ow=-Pshrink correspond to zero slope, i.e. extrema of 
the free energy. Crossings of P gr0 w and -P s hrmk, then, can be used to determine 
the locations of the extrema, including the saddle point. 

Here we present results for Projective Dynamics applied to micromagnetic 
simulations of magnetization reversal in a chain of 17 spins, Si, with the chain 
aligned along the z-axis [8] . Using an extended Heisenberg model, the internal 
energy of the system is given by 
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where J is the exchange energy, B is the strength of the external field oriented 
parallel to —z, and D is the strength of the dipole-dipole interactions in energy 
units, see Eqs. (17) and (18) of Ref. [9]. We choose parameters consistent with 
previous studies of iron nanopillars modeled as a one-dimensional chain of 
spins [7,8,9,10]: J=1.6 x 1(T 12 erg and D=A.l x 1(T 12 erg. The dynamics 
consist of each spin precessing around a local field, i.e. the Landau-Lifshitz- 
Gilbert (LLG) equation [11,12] 
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where the scaled electron gyromagnetic ratio 70 = 9.26 x 10 21 Hz/erg for 

— * 

this system, and the local field is given by the functional derivative Hi = 
—SE/SSi. The phenomenological damping parameter a = 0.1 was chosen to 
give underdamped dynamics, and an Euler integration scheme was used [13]. 

Dipole-dipole interactions make head-to-tail alignment of the spins along the 
chain favorable and provide a strong uniaxial anisotropy for this system. For 
each simulated reversal, the nanopillar was allowed to come to equilibrium in 
a field B —1.9 x 10~ 12 erg (1000 Oe). Then the magnitude of B was quickly 
decreased to a value of — B to create a metastable configuration of the mag- 
netization. For each integration step, data for Pg roW ) -^shrink, and the energy 
were binned on M z . Results for this system have been presented previously 
[7]. Here we present improved results obtained from a systematic choice in bin 
width and from including lower temperatures. 

Projective dynamics is most easily analyzed for discrete systems, such as the 
Ising model with spin-flip dynamics. Then M z changes by —S, 0, or +S at 
each step, and the stochastic dynamics are easy to analyze in terms of jumps 
along a well-defined chain of states. When M z is a continuous variable, binning 
is used to create the chain of states. However, there are competing limitations 
that determine the bin size. One is the assumption of complete mixing within 
a bin, which is required for the transitions between bins to be a Markov pro- 
cess. If the bins are too large, then complete mixing within the bins will not be 
present and there will be memory effects. Another limiting effect is gathering 
enough statistics in each bin to measure the probabilities P gr0 w and P s hrink- If 
the bins are too small, there will be inadequate statistics related to these mea- 
sured probabilities. Furthermore, we would like to have the projected Markov 
matrix be tridiagonal, which allows calculation of residence times using Eq.(3), 
introduced below, rather than inverting a large matrix. Tridiagonality of the 
Markov matrix requires that jumps from anywhere in the bin can go only into 
the bin with the next largest or smallest ^-component of the magnetization. 

The results presented in this work were obtained in the following manner. Data 
were collected for on the order of 10 2 switches for each temperature with a bin 
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Fig. 1. Probability of shrinking, P s hrink) (squares) and growing, P gr0 w; (circles) for 
T=5 K. The solid curves are fifth-order polynomial fits to the data. The locations 
of the metastable minimum and saddle point, as indicated by the arrows, are deter- 
mined from the intersections of the polynomials. 

of width 6=5 x 1CT 5 . The data were then analyzed for bin widths of 6, 26, 36, 
etc. until jumps only occurred between neighboring bins. Using this method, 
the optimal bin width was found to be 56 at T=3 K and 256 at T=50 K. In 
contrast, in Ref. [7] bins of width 1006 where used. Both probabilities, P gr0 w 
(circles) and P sn rink (squares), for T=5 K, are shown in Fig. 1. Also shown are 
two fifth-order polynomials fit to these probabilities. Throughout this work, 
fifth-order polynomial fits have been used to determine the intersection of 
-fgrow and P s hrink, and we find that the order of polynomial used affects the 
estimated location by less than 1% for the saddle point. The crossings for 
T=5 K are labeled with arrows in Fig. 1, with the crossing near M 2 =0.995 
the metastable minimum and that near M z =0.95 the saddle point. 

The values of M z for the metastable minimum (circles) and the saddle point 
(diamonds) for all temperatures T are shown in Fig. 2. Similar results have 
been shown previously [6,7], but the results shown here are quantitatively 
different because of the smaller bins used. The changes are as large as 10% 
and are strongest for the saddle point. The linear dependence of the M z of the 
saddle point on T noted previously [6,7] is present in the new data as well, 
but only for T>5 K. The solid line is a least-squares fit, and its intercept is 
M 2 =0.958. The deviation from linearity at low temperature is expected to the 
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Fig. 2. M z of metastable minimum (circles), internal energy maximum (squares), 
and saddle point (diamonds) vs temperature. The separation between the saddle 
point and the internal-energy maximum implies large entropic effects in magneti- 
zation reversal. The solid line is a least-squares fit to the saddle-point data. The 
metastable minimum is identically unity at T=0 K. 

extent that extrapolation of the linear fit to T=0 K is not consistent with the 
internal-energy maximum, described below. This deviation was not seen in 
previous studies since such low temperatures are considered here for the first 
time. 

The free energy, F=E — TS, has contributions from the internal energy and 
the entropy, S. When the temperature is exactly zero, there is no entropic 
contribution to the free energy, and the saddle point corresponds to the max- 
imum in the internal energy. The results for the maximum of the average 
internal energy, (E), binned on M z are shown (squares) in Fig. 2. The M z for 
the saddle point and the internal-energy maximum converge as the tempera- 
ture decreases, as expected. The separation between M z at the internal-energy 
maximum and the saddle point at fixed T is a measure of the importance of 
entropic contributions to the magnetization reversal process. The separation 
in Fig. 2 indicates that entropy is quite important in these model nanopillars. 

The internal-energy "barrier" , the difference between (E) at the saddle point 
and the metastable minimum, is shown vs temperature T in Fig. 3 as squares. 
Note that for all but the lowest temperatures the difference is negative, indi- 
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Fig. 3. The free-energy barrier (circles) and internal-energy barrier (squares) vs tem- 
perature. The difference between the two barriers is due to entropic contributions 
and vanished as T^0. The difference between the maximum internal energy and 
the internal energy at the metastable minimum (diamonds) is also shown. 

eating that the (E) for the saddle point is lower than that of the metastable 
minimum. It is indeed possible for a system to be metastable in this situation, 
because it is a free-energy barrier that is actually required. The free-energy bar- 
rier can be estimated using Projective Dynamics from the time spent in each 
bin i, h(i), which is proportional to the Boltzmann factor exp(— F(i)/kBT). 
For a tridiagonal Markov matrix, the residence time h(i) can be found from 
the growing and shrinking probabilities by [4,5,3] 

h{i) = [1 + P shrink (i - l)h(i - 1)] /PgrowW , (3) 



with h(l) = l/P grow (l), where we take the first bin to be the stable state 
and iterate Eq. (3) towards the metastable state. (Note that the cited ref- 
erences iterate from the metastable to stable state.) free-energy barrier is 
AF=/c s Tln(/i(v in )//i(i sadd i e )), where i min is the bin containing the metastable 
minimum and i sa ddic contains the saddle point. The results for the free-energy 
barrier are shown (circles) in Fig. 3. The free-energy barrier is always positive 
and approaches the internal-energy barrier as T— >0, as expected. However, it 
is clear from the internal-energy barrier being negative for most temperatures 
and from the temperature dependence of AF, that the entropy contributes 
significantly to the metastability. This indicates that the configuration of {Si} 
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is relatively restricted at the saddle-point, and that this restriction is more 
important than the combined dipole-dipole and exchange energies at all but 
the smallest temperatures. The difference between the maximum (E) and the 
value at the metastable minimum is also shown in Fig. 3 as diamonds. This 
difference might easily be mistaken for the barrier to magnetization rever- 
sal, but here it is much less sensitive to the temperature than the measured 
free-energy barrier. 

In this paper, the usefulness of Projective Dynamics for locating the metastable 
minimum and saddle point in magnetization reversal has been demonstrated. 
Once these points are known, it is also possible to directly find the free-energy 
barrier of the reversal process, also using Projective Dynamics, or the internal- 
energy barrier, using some other technique to find the internal energy. One 
of the truly exciting aspects of Projective Dynamics is that it is not limited 
to analyzing simulations. Since only data on the time- dependence of the slow 
variable is needed, the analysis can also be applied to experimental data. 
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